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C/^ ; Abstract 

\l ' We clarify the procedure for expressing the Friedmann equation in terms of directly measurable 

^^ ■ cosmological scalars constructed out of higher derivatives of the scale factor. We carry out this 

procedure for pure dust, Chaplygin gas and generalised Chaplygin gas energy-momentum tensors. 
O , . In each case it leads to a constraint on the scalars thus giving rise to a test of General Relativity. 

We also discuss a formulation of the Friedmann equation as unparametrised geodesic motion and 



^ . its connection with the Lagrangian treatment of perfect fluids coupled to gravity. 



CN ■ 1 Introduction 

>, 

It is a striking and slightly puzzling fact that almost all current cosmological observations can be 

^ . summarised by the simple statement [H [21 [12] 

O 

The jerk of the universe equals one. 



The aim of this paper is to uncover the mathematical reason behind this fact and explore some of its 



o 

Q , consequences, and generalisations. 

Our present universe appears to be well described by the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) metric 

g = -dt^ + a{tfh, (1.1) 



where the function a = a{t) is the scale factor, /i is a metric on H"^, W^ or S^ with constant curvature 
/c = — 1, or 1 respectively and the speed of light has been set to 1. 
One may define four cosmological scalars by 

Ida fda\~'^<f'a ^ 2(da\-^d^a ■^/da\-'^d'^a 



These scalars are known as Hubble, deceleration, jerk and snap respectiveljo. The last three of these 
are dimensionless and the Hubble has the dimension of inverse of time. The scalars can in principle 
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^ The analogous expressions involving 5th and 6th derivatives of a are known as crackle and pop. This terminology 
goes back to a 1932 advertisement of Kellogg's Rice Crispies which 'merrily snap, crackle, and pop in a bowl of milk'. 



be measured using red-shift data and the Hubble law [3 El [H [121 [II [21 [13]. They arise, for example, 
as the first four terms in the Taylor expansion of the red-shift 1 + z = a{to)/a{te) observed at time 
t = to as a function of time of emission t = te- It is important to stress that the Einstein equations 
have played no role so far. The form of the FLRW metric is derived from symmetry consideration 
which requires no knowledge of the dynamical variable o(f). The Hubble law is an observational fact 
also independent of the field equations. 

In this paper we shall regard the Einstein equations, in the form of Friedmann equation for a = a{t), 
as one algebraic constraint between the scalars (II. 2p . This links the measurement of the cosmological 
scalars to a test of General Relativity, or any of its modifications (which would lead to different 
constraints). If one assumes that Einstein equations hold then measuring the cosmological scalars 
could determine the the equation of state in relating the energy density and the momentum in the 
perfect fluid energy momentum tensor. We aim to clarify this procedure (which appears to have been 
initiated in [7j and developed by fSj and is closely related to the Statefinder approach of [l2]) and link 
it to some recent work [3] on the metrisability of projective structures. For example it turns out that 
the constant jerk condition Q = 1 is equivalent to Einstein equations with k = and the dust stress 
tensor. Using the scaling and translational symmetries of the ODE Q = 1 allows to put the general 
solution in the form 

a(t) = sinhi (y^i), 

where A is the cosmological constant. Another motivation comes from quantum cosmology. In the 
Hamiltonian treatment of general relativity one regards the Einstein equation as a dynamical system 
equivalent to a forced geodesic motion on an infinite dimensional space M of three-dimensional 
Riemannian metrics on initial data surface S. The metric on Ai is the DeWitt metric 



0,(^9, Sg) = J VaiM^g^) - Tr{6gf)d^x 



(where 6g £ TgM and the traces are taken with respect to g) and the potential is given by the scalar 
curvature of the Riemannian metric on E. See e.g. [5] where this is discussed. In the context of 
quantum cosmology this analysis needs to be modified in the presence of matter and cosmological 
constant. In Section [3| we shall use a different point of view to relate the Friedmann equation a 
geodesic motion on two-dimensional surface of revolution. In Section [3| we shall make contact with 
the Lagrangian treatment of perfect fluids coupled to gravity. The problem of finding a metric whose 
unparametrised geodesies coincide with the integral curves of a second order ODE is summarised in 
the Appendix. 

2 Cosmological scalars and the Friedmann equation 

Consider the Einstein equations with cosmological constant A and matter described by the perfect 
fluid energy momentum tensor 

T^,u = {p + p)U^U^+pg^^ (2.3) 

where f7^ = V^t, p = p{t) describes the matter density and p = p{t) is the pressure. The energy- 
momentum conservation yields 

dp , .da 

and the Einstein equations reduce to the Friedmann equation 

da\2 8-kG n a n ,„ _,, 

_) +t=—pa-^-a^ (2.4) 



Pure dust. Consider the pressure-free matter p = 0. The energy conservation gives pa^ = M 
where M is a constant. Let us now consider a system of three equations consisting of ()2.4p and its first 
two time derivatives. We regard this as a system of algebraic equations for the constants {k, A, GM) 
which can therefore be expressed as functions of (a, a, a, 'a). Take the third derivative of (j2.4[) and 
substitute the expressions for {k,A,GM). The resulting equatioro does not contain any parameters 
and can be expressed in terms of the cosmological scalars p.2p as 

X + 2{q + Q) + qQ = 0. (2.5) 

This fourth order ODE is equivalent to the Friedmann equation and has an advantage that it appears 
as a constraint on directly measurable quantities. Thus it provides the test of the model alluded to in 
the Introduction. If only two constants (A, GM) are eliminated between (j2.4p and its first derivative 
then the second derivative of (j2.4p yields 



k = a'^H'^{Q - 1) 

where k is regarded as a parameter. This is the formula obtained in [7]. In particular if A; = this 
relation reduces to a third order ODE 

Q = l. 

This constant jerk condition is consistent with recent redshift analysis which provides the lower bound 
for the jerk implying that Q > 0. Compare a related discussion in |13j . 

Chaplygin Gas. This exotic form of matter is given by the perfect fluid energy momentum tensor 

with the equation of state 

A 

P = , 

P 

where ^ is a positive constant. In the cosmological context this model was proposed in [9]. The 

negative pressure allows to describe a transition from a decelerated universe to cosmic acceleration. 

The energy conservation equation gives 



p = a/^ + |- 



For small a{t) this reduces to the dust cosmology p = \/Ba ^, and for large a one gets the de Sitter 
Universe p = vA,p = —yA. In between these two regimes one can use the expanded expression 

p = Va+ ^ -« 



2^/Z 



Thus y A plays the role of a cosmological constant. We insert this to the Friedmann equation with 
A = and follow the procedure of eliminating the constants by differentiation. This leads to an 
approximate constraint 

X + b{Q + q)+qQ = Q (2.6) 



^This idea of eliminating parameters and reinterpreting them as constants of integration has a long history. It goes 
back at least to Halphen [Bj who obtained a fifth order ODE characterising conies in RP^. In the inhomogeneous 
coordinates (a, t) the five parameter family of conies is 

2 2 

a — c\ t + C2 ai + C3 a + C4 i + cs . 

Eliminating the parameters (ci,...,C5) between this equation and its fourth derivatives and substituting in the fifth 
derivative yields the ODE 



which is vahd when a ^B l2\fA is small compared to the unity. 

Generalised Chaplygin Gas. This is the generalisation of the previous case, where 

P = -4' 0<a<l. 

P 

Introducing an additional constant a allows a description of a universe evolving from the pure dust 
matter to a cosmological constant with an intermediate epoch with non-zero cosmological constant 
and perfect fluid matter satisfying the equation of state p = ap. 
The energy conservation gives 

1 

p={A + i?a-3("+i))^ ^A^ + ^i:li4a-3("+i) + . . . . 

a + 1 A 

Inserting the expanded formula into the Friedmann equation with A = leads to an expression with 
three arbitrary constants. Our procedure gives the constraint 

X + {3a + 2){Q + q)+qQ = 

in agreement with ()2.5p and ()2.6p when a = and a = 1 respectively. We may however want to 

eliminate a from this equation by another differentiation. This requires an additional cosmological 

scalar crackle 

4 /(ia\-5(i^a 



\dtJ dt^ 
A rather complicated (MAPLE) calculation leads to a relatively simple constraint 

-2qX - 2Qq^ -Yq- 2XQ - 3Xq^ - Q^q -YQ + X^ - 3q^Q - qXQ + Q^ - 2Q^i = 

which should hold if Einstein equations with the generalised Chaplygin gas energy momentum holds. 
This constraint is again approximate and is valid only in the regime where the higher order terms in 
the expansion of p can be dropped. 

A different constraint involving Y appeared in the work of Hut [8] who considered a mixture of 
perfect dust and radiation in the energy momentum tensor. 

3 Friedmann equations as geodesic motion 

In this section we shall relate the Friedmann equation to a geodesic motion on a surface of revolution 
with metric given explicitly by (j3.8p . 

For simplicity we consider the pure dust case p = 0. Eliminating the constant GM between the 
Friedmann equation and its time derivative yields a second order ODE 



^^ = UAa-^A--(-S'. 



f3 7) 
dt"^ 2V""' aJ 2a\dtJ ' ^ ' ' 

Can we think of the integral curves of this equation as unparametrised geodesies of some metric on 
a surface with local coordinates (a,i)? This question can be asked about any second order ODE of 
the form a = F{t,a,d). The necessary condition is that F is at most cubic in the first derivatives d. 
This condition (which was known to R. Liouville [10]) is obviously satisfied by ()3.7p . Therefore the 
integral curves of (j3.7p are geodesies of some projective structure (an equivalence class of torsion-free 
connections sharing the same unparametrised geodesies). In [10] Liouville also began the study of 
sufficient conditions for a projective structure to come from a metric, but these were derived only 



recently [3] . They come down to vanishing of three weighted invariants of differential order at most 8 
in connection defining the projective structure. We have checked that these invariants vanish of (13. 7|] 
and so there exist an underlying metric. Following a linear algorithm lied down in [lOj and summarised 
in the Appendix one can construct the metric explicitly. The answer is 



ka-^a^ + c {ka- ^a^ + cy'' 

where c is a constant equal to a negative multiple of the Newton's constant. To verify this, write 
the geodesic equations for (jS.Sp and eliminate the affine parameter between the two equations thus 
expressing a as a function of t. This will lead back to (j3.7p . The scalar curvature of the metric ()3.8p is 

2kAa^ + 9Aa^c + 6k^a + 3kc 

Thus the metric has constant curvature (and therefore it is projectively flat) if and only if /c = 0. 
In this case there exist a coordinate transformation (t, a) -^ {t(t,a),d{t,a)) such that, in the new 
coordinates the Friedmann equations are 

d'a 

IP-"- 



The cubic denominators in (j3.8p have real zeroes, and so the metric may appear singular. In fact these 
zeroes correspond to coordinate singularities, and are not present in curvature scalars. For example 
when k = the curves a^ = 3c/ A are only apparent singularities as the metric has constant curvature. 

4 Lagrangian description of perfect fluids 

In this section we shall make contact with a standard Lagrangian treatment of perfect fluids coupled 
to gravity [11] and then specialise to the FLRW case. We shall then use this formalism to derive the 
geodesic formulation analogous to (13. 8p when radiation and other forms of matter are present in the 
energy-momentum tensor. 

We shall restrict our general discussion to the case of an irrotational (i. e. vorticity-free) perfect 
fluid. We introduce a potential ^ and let 



We take as Lagrangian for ip 
The energy momentum tensor is 



L = L{Z) . 



Comparing this with a perfect fluid (j2.3p with U^ = (SJ ^'il))/Z we may identify 

p = 2ZLz — L p = L . 
Adding a constant to the Lagrangian L(Z) introduces the cosmological term in the FLRW models: 
L^L + 7 p^p-j, V^P + 1, Tf,^ — >Tf,^ + 7gii,y 
• If p = ap then L = Z ^a . 



• For a Born-Infeld scalar field 

L = -Vl-^ + l 

we get 

P 

Omitting the one, which amounts to changing the zero of the energy scale, or cancelling a 
cosmological term, gives the Chaplygin gas 

1 

P = — 

P 

mentioned earlier. 

The shift symmetry ^ — > ^ + constant gives rise to a conserved current which, in this perfect fluid 
context, may be interpreted as the entropy current. 

As an example, consider radiation L = Z^. The equations of motion are 

V^((VV)'V^^) =0, (4.9) 

or, as long as V^tp is time-like 

Two simple solutions in flat space-time are 

• tp = fit — z), where z is one of the spatial coordinates and / is an arbitrary function. This gives 
a trivial solution with vanishing energy momentum tensor. 

• ip = t. This represents a uniform fluid at rest. Linear perturbations about this, or indeed any, 
solution are easily seen to travel with respect to it with speed -h=. 

Interestingly the action, and hence the equations of motion, are invariant under Weyl conformal 
rescalings 

V'^V', g^il^ ^^'^{x)gf^u■ 
ln other words the equations (j4.9p are conformally invariant, just as is Yang-Mills theory in four 
space-time dimensions. 

4.1 Coupling to FLRW models 

To couple to gravity we would consider the Lagrangian ( up to a boundary term) 



+ L(Z) V=5d2;. (4.10) 



.IGnG 

To obtain the full consequences of the Einstein equations we substitute the ansatz ip = 'ip^t) and 

g = -N'^df+a{tfh, (4.11) 

where A'^ = N{t) is an arbitrary lapse function into the action ()4.10p . We vary with respect to ip^a 
and A^. The first two Euler-Lagrange equations are second order in time. The third variation gives 



a constraint on the initial values ■0) ^0 and a, a which is consistently evolved by the two second order 
equations. It is convenient to make the choice A'^ = 1 after variation. It then follows that the constraint 
is just the vanishing of the Hamiltonian 7i obtained from the Lagrangian 

C = —L{Z) -{a{hf-ka), where Z = 0^ (4^2) 

4 SttG 

by a standard Legendre transform. One may check that the resulting system of ordinary differential 
equations is the same as the system one would obtain by substitution of the ansatz (j4.1ip into the full 
Einstein equations. 

Clearly, for a general equation of state p = p{p), the kinetic term in the Lagrangian for ^ will not 
be quadratic in ip and hence the standard Jacobi procedure will not lead to a metric. The exceptional 
case is of course 'stiff matter', p = p which is well known to be represented by a standard massless 
scalar field. We shall instead use the procedure introduced in Section [3] to construct a metric whose 
geodesies coincide with unparametrised extremals of the Euler-Lagrange equations of (j4.12p . 

The conservation of momentum a?il)L'{Z) =const can be used to find ip = ipla). It is convenient 
to introduce a function f = f{a) such that 

^ = 27rGa^L{Z{a)) + k. 
da 

The equations of motion resulting from ()4.12|) become 

^ = -J-^-^f^f (4 13) 

dt'^ 2a da 2a\dt) ^ ' 

Examining the obstructions of [3] shows that the projective structure defined by (J4.13P is metrisable 
for any f{a) (and therefore for any choice of the Lagrangian L). Following the algorithm of [3] leads 
to the expression for the metric 

dt^ a da^ 

Unparametrised geodesies of this metric are integral curves of (j4.13p . For example \i p = ap and the 

l + CK 

cosmological constant is allowed then L{Z) = Z^^ + 7 and 



^ = 777Y + T732- (4-14) 



2-KaG f b \i+"_^3„_i 27rGj 3 



where b and c are some constants. The pure dust case considered in Section [3] is recovered in the limit 
a — > 0. In this case / is cubic in a and the metric (I4.14|) becomes (13. Sh if 7 = — A/(27rG). 
For radiation (a = 1/3) the Legendre transform gives the Hamiltonian and the momenta 

n = ^ii^)' - ^{a{df + ka) , 
p^ = a^{ipf, 

where H = d/o. Note that the conserved momentum p^ is essentially the conserved entropy of the 
radiation (i.e. number of photons). The phase ip is ignorable both in the technical and the vernacular 
sense since it disappears from the Hamiltonian. 



5 Summary 

We have studied the Priedmann equation from the point of view of geometry of ODEs. 

In Section 2 we set up an algorithm replacing the 2nd order Friedmann equation which involves a 
number of parameters (this number depends on a choice of equation of state) by a constraint on directly 
observable cosmological scalars. The validity of resulting constraints can in principle be verified by 
observations, thus providing the experimental test of General Relativity as well as railing out some 
equations of state. For example the Priedmann equation with k = and the dust energy-momentum 
tensor is equivalent to the constant jerk condition Q = 1. 

The constraints discussed in Section 2 give rise to curves in the (a, t) plane, one curve through 
each point in each direction, and in Sections 3 and 4 we have demonstrated that these curves are 
unparametrised geodesies of some two-dimensional (pseudo)-Riemannian metric. This is an alterna- 
tive to the more usual DeWitt approach and may play a role in projective quantisation of Cosmology 
which does not refer to a specific (time) parametrisation of classical trajectories. 

Appendix 

Recall that a projective structure on an open set [7 C K^ is an equivalence class of torsion free 
connections on TU. Two connections are projectively equivalent if they share the same unparametrised 
geodesies. Let (a, t) be local coordinates on U. Eliminating a parameter s between the two geodesic 
equations x* + T\x^x^ = ux* where x*(s) = (t(s),a(s)) yields a second order ODE 

^ = ^3(M)(f)V^.(M)(|)VAi(M)(f)+^o(M) (Al) 

where the functions Aa are given in terms of the coefficients of a connection in a given projective class 

^o = -r?i, Ai = r}i-2r?2, A2 = 2t\^-tI^, A3 = r^2- (A2) 



Conversely, any ODE of the form (jAip gives rise to a projective structure. 

A projective structure is called metrisable if the integral curves of the ODE (jAip are unparametrised 
geodesies of a metric 

g = E{t, a)df + 2F{t, a)dtda + G{t, a)da^ . (A3) 

In this case the four functions A^ are given in terms of (-E, F, G) and their first derivatives. The 
corresponding expressions can be derived by calculating the Levi-Civita connection of g and reading 
off the AnS from (IA2I). The substitution 



linearises these expressions thus leading to a characterisation of the metricity condition as consistency 
conditions for an overdetermined system of linear PDEs 

Theorem 1 (R. Liouville 1889 [lOj). A projective structure corresponding to the second order ODE 
(jAip is metrisable on a neighbourhood of a point p £ U iff there exists functions tpi,'4'2,'^3 defined on 
a neighbourhood of p such that 



does not vanish at p and such that the equations 

di^s di;2 . , 4 ^ 2 , , . , 

J^ + 2^ = 2A,.p,--A,^, + -A2i.2 (Al) 

/loW on i/ie domain of definition. 

This system is overdetermined, as there are more equations than unknowns. In [3] the consistency 
conditions for this system where found in terms of point invariants of the associated second order 
ODE (jAip . The details of this construction and the invariants themselves are rather complicated and 
we refer the reader to [3]. The invariants obstructing metricity vanish for the projective structures 
()3.7p and ()4.13p arising in our paper, and the corresponding metrics are found by solving the linear 
system ([Mj) . 
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